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A complete set of asymptotic three particle light cone distribution amplitudes of 
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projectors technique which allows to handle kinematical corrections due to a finite 
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production and 7p7r form factor are discussed. Our results may imply a breakdown 
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submitted to Phys. Lett. B 



1 Introduction 



Light cone distribution amplitudes (DA's) of hadrons are the key ingredient of the 
modern QCD approach to hard 0] and soft exclusive , diffractive and deep 
inelastic scattering [Q processes. They provide the universal non-perturbative input 
in physical amplitudes where the only dependence on a particular process enters as a 
normalization point of a corresponding DA. Respectively, they have to be addressed 
using non-perturbative methods. The approach of the QCD sum rules pioneered for 
the study of light cone DA's by Chernyak and Zhitnitsky more than a decade ago |^ 
hitherto remains the most popular and reliable among methods intended to handle 
this problem. 



While a good deal of information was obtained on leading twist DA's for different 
hadrons, relatively a little has been known about nonleading twist DA's except for 
the case of pion. A systematic approach to the study of higher twist DA's was 
suggested in where a complete set of light cone DA's of twist 3 and 4 for the 
pion has been constructed and studied in detail. The analysis of Ref . ^ was heavily 
based on group-theoretical methods exploiting the conformal symmetry which holds 
in QCD at the one-loop level l]^? 0- The conformal invariance has proved useful 



for the study of leading twist DA's [|10| where it allows to diagonalize the mixing 
matrix at one loop. The use of the conformal symmetry becomes even more crucial 
in the analysis of nonleading twist DA's where methods based on direct calculations 
of suitable correlation functions turn out either too complicated or erroneous p. 



The aim of this letter is to extent the conformal group methods of Ref.^ to 
a study of nonleading twist three particle light cone DA's for a massive vector 
meson built out of massless current quarks Q Fortunately, it requires only a slight 
modification of the technique developed in to handle kinematical corrections due 
to a non-vanishing hadron mass. This method yields quite a different classification 
and new results for twist 4 DA's in comparison to the analysis of Ref.[^. Our 
interest in this problem is mainly due to its utility for the diffractive vector meson 
electroproduction at HERA |T2| (see |]T3| for updated review of the experimental and 



theoretical state of affair there). Recently the perturbative QCD approach of Ref. 
to the diffractive electroproduction of a longitudinally polarized vector meson was 
reformulated directly in the coordinate space in the language of Operator Product 
Expansion (OPE) The method of Ref. easily incorporates an asymmetry 



of the gluon distribution in the nucleon which is important for the whole problem 
| |15| |. Moreover, higher twist DA's appear quite naturally in this approach. Though 
this paper deals with a transversely polarized vector meson for which the method 



of [|14| must be somewhat modified, as a by-product of our study we show that 
there are no three-particle DA's in the diffractive amplitude of Ref. p, Q ( see 
the end of Sect. 3). Among other potential applications one can mention the light 



"See Ref.||ll| for a recent re-analysis of leading twist p-meson DA's. 
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cone QCD sum rule approach to semileptonic and radiative B-decays |]T^ and 
(non-)factorization in non-leptonic B-decays ||18[|. Besides these practical issues, our 
results may have some conceptual importance as they seem to imply a breakdown 
of OPE for some exclusive processes at the nonleading twist level. The letter is 
organized as follows. In Sect. 2 we review the conformal group technique of Ref.p. 
Its modification for the case of interest is presented in Sect. 3 where we construct a 
complete set of asymptotic three particle DA's for a transversely polarized massive 
vector meson. Sect. 4 contains a comparison with previous analyses and discussion 
of possible consequences of our results. 



Conformal symmetry on the light cone and 
asymptotic distribution amplitudes 



The aim of this section is to remind the algebra of the so-called collinear conformal 
group obtained by the reduction of the full conformal group to fields varying on the 
light cone = and establish a projection operator technique and general formulae 
for asymptotic distribution amplitudes with arbitrary number of constituents. 

The algebra of the conformal group 0(4,2) is obtained by adding the generators 
D of dilatations and of special conformal transformations to the generators 
and M^y of the Poincare group. They form the algebra 

[D,M^,] = , [K^,K]=0 (1) 
[K^,P,] = -2i{g^,D + M^,) 
[Ka, M^^] = i{gaf,K^ - gauKf,) 

in addition to the commutation relations of the Poincare algebra. These generators 
act on a field $(a;) with the spin s and canonical dimension / as follows 

[P^,$(x)] = -td^^x) 

[D,^x)] = -lixada + lMx) (2) 
[K^,(^{x)] = -i{2Xf,Xada- x'^d^ + 2Xf,l -2ix^J:^^)(^{x) 

where T,^i, is the spin generator 

^fiu"^ = 2^^;^^ ' ^tiuGap = i{S^aGui3 - SuaG^is) - {a ^ P) (3) 

For an analysis of the light cone hadron distribution amplitudes we only need to 
consider the action of the conformal generators (|1]) on a field $(-ux) varying on the 
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light cone x = 0. It is convenient to simplify the algebra by introducing two light like 
vectors x^, x^ with = = , (xx) = 1 and defining projections on these vectors 
7. = 7^a;^ , 7^= = 'jfj.Xfj^ etc. Then the only nontrivial transformations on the field 
are given by the components P. , D , M*. , K^,. These generators form the 
algebra of the so-called coUinear conformal subgroup 5*0(2, 1) ~ SU{1, 1) ~ SL2{R) 
of projective transformations on the line. Its algebra reads f\ 

[J+,J_] = -2J3 , [J3,J±]=±J± 

= , [E, J,] = (4) 



where we defined the operators Ji as hnear combinations of the conformal operators 

J+ = Ji+ iJ2 = ' J- = Ji — iJ2 = y^-^* 

^3 = ^P + M,.), E='-{D-M,.) (5) 

(/ = t/g — 'J^ — 'J2 ~ ^3 — — ^+ 

Here two observations are in order. First, the light like vectors , x^ introduced 
above are quite arbitrary at this stage. However, they have to be built up from 
physical vectors characteristic to a meson in order to get informative constraints 
for distribution amplitudes. Ways to do this will be described in the next section. 
Second, we note that the use of the coUinear conformal group requires fixing a 
projection of the Lorentz spin of the field $(mx) onto the line x^ 

S^.<I>(mx) = ^^^(x) (6) 

In this case the field $(mx) is an eigenstate of the Casimir operator : 



J , ^(ux) 



j(j - iMux) (7) 



where j = l/2(/ + s) is the conformal spin of the field $(mx). The fixed spin 
components of the quark \E'(-ux) and gluon G^i,{ux) fields can be readily constructed 
as follows. The quark field contains two components with the spin projections 
s = ±1/2 which can be selected by the projection operators 7. and 7=,, : 

S,.7.^' = ^a*.7.* = +z^7-^ 

S*.7*^ = ^fT*.7,* = -z^7*^ (8) 

For the gluon field the relevant components with the spin projections s = 0, ±1 are 
the following 

Tn:.G.± = -\-iG.± , 12^.G^:± = —iG^± , S^.G^. = T,^.G±± = (9) 



^Note the "wrong" sign of the commutator [J+, J_] in eomparison to the rotation group S0(3). 
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In the last equation the label _L stands for projections on the third vector e^^^ which 
is required to satisfy (xe^-^^) = 0. 

In an analysis of hadron DA's we deal with matrix elements of the form 

(0|$i(Mix)...$fcKx)|%)) = J Z}ae-^''^("^"^+-+"^"'=V(ai,---,afc) (10) 

where Da = dai ■ ■ ■ ak^iJ^i c^i ~ 1)- Here /i is a hadron of interest and 0(ai, . . . , ak) 
stands for a corresponding multiparticle DA. The fields ^i{uix) are arranged along 
the light cone = 0. The light cone conformal symmetry suggests the following 
way for the analysis of the DA 0(ai, . . . ,afc)- First, we choose fixed spin compo- 
nents of the fields ^i{uix) according to Eqs.@,(^). Next, we note that an expan- 
sion of the tensor product ^i{uix) . . . ^k{ukx) into irreducible representations of 
the coUinear conformal group induces a respective decomposition for the function 
. . . , ak) defined on the simplex ai + . . . + ak = 1. Therefore the evaluation 
of the DA (f){ai, . . . , ak) amounts to a calculation of corresponding Clebsch-Gordan 
coefficients in an expansions of a tensor product into irreducible representations with 
a fixed value of the total conformal spin. An asymptotic DA is then defined as a 
representation with the lowest conformal spin which is jmin = Ji + • • -jk- A general 
formula for the asymptotic form of arbitrary multiparticle DA has been found in 
Ref. exploiting the fact that the collinear conformal group (H) is identical to the 
Lorentz group in (2+1) dimensions. It reads 

= Y^^— ^ f(2j^"^ ' J Da(l)as = l (H) 

Irreducible representations with higher spins are obtained by the multiplications of 
(|ll|) by mutually orthogonal polynomials of dimensions j — J2i ji ^-nd correspond to 
conformal corrections to the asymptotic multiparticle DA (j){ai, . . . , ak) which serves 
as the weight function at the space of conformal polynomials. These conformal 
corrections will not be discussed in the present paper. 

Going over to interacting fields we insert the gauge factors [ux, vx] = 
Pexp[ig dtXf^Af^itx)] in between the field operators in Eq.(p!ol) which is equivalent 
to the substitution of ordinary derivatives by covariant ones in the language of local 
conformal operators. Then contributions of different conformal spins in DA do not 
mix under renormalization to the one- loop accuracy [0. This property is sufficient 
to construct an expansion of the leading twist DA over a series of polynomials with 
multiplicatively renormalizable coefficients, as in this case all representations are 
nondegenerate. Moreover, anomalous dimensions are ordered with the conformal 
spin for the leading twist. For higher twist DA's irreducible representations with 
sufficiently high spins are degenerate and conformal operators do not completely 
diagonalize the mixing matrix which only becomes "conformal block-diagonal" . Still, 
anomalous dimensions are typically lowest for operators of minimal conformal spin. 
These operators are therefore leading ones as — ^ oo and thus in this case the 
definition (0) is justified. 



4 



3 Distribution amplitudes of a massive vector me- 
son 



In this section we explain how the general analysis of the proceeding section works in 
the study of three-particle quark- ant iquark-gluon DA's of twist 3 and 4 of a massive 
vector meson. We consider a task of taking into account kinematical corrections 
due to a nonvanishing meson mass while the constituents of a meson are assumed 
to be massless. Thus, our results apply directly to the p— , cj— , (p— mesons (SU(3) 
limit is implied) and also can be used for the charmonium in situations where one 
can neglect the c-quark mass. For definiteness, we consider the case of the charged 
p— meson, while modifications to be done for other mesons are trivial. 

We start with a set of three particle DA's of twist 3 and 4 defined by the matrix 
elements of the following non-local light cone string operators [| 

{Q\d{-x)-iai9G^u{vx)u{x)\p{q)) = Qais^^^u - ei^kf.) J D{aiqx)^^^\a^) 

2 

+£a\qt.Xu - QuX^,)— f D(aigx)$(2)(a.) (12) 
qx J 



{0\d{-x)jal5gGi^u{vx)u{x)\p{q)) = qa{e^^hu ~ j D{aiqx)¥^\ai) 

2 

+£l^Kqf.x^-quXf,)— [ D{aiqx)¥^\ai) (13) 
qx J 

+ {ei'^Sau-ei''>6^,)m' J D{a,qx)¥^\a,) 

Here 

J D{aiqx) = J daida2das6{l - ai - - a3)e-^^^("^-°^+''"^) (14) 

and stands for the meson mass. One can check that the DA's corespond 
to operators of twist 3, while the DA's '^^^^ contain contributions of 

operators of both twist 3 and 4. This can be seen using the operator identity M 



tw-3 1 + d 

2 dxo, 



i^{-x)'ialb9G^.{vx)i){x) = I du — - — ^^il){-ux)'^.-i^gGf,.{uvx)i){ux) 



/■I X — - 
+ du — - — —ip{-ux)-f.-f5gGa.{uvx)ilj{ux) (15) 

Multiplying (^) by Xa, integrating by parts and taking matrix elements in accord 
with (|13|), one arrives at identities for the DA's <l>(^\ \E'(^\ We further note that the p- 
meson in Eqs. (|T2|) ,(|l3D can only be transversely polarized. Indeed, the total helicity 

^We assume the Fock-Schwinger gauge Xf_iA^{x) — 0. In a general case path-ordered gauge 
factors between the fields in (|T2,0) are implied, as required by the gauge invariance. 
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of the quark-antiquark pair is zero for the vector or axial channel. Therefore, in 
order to have a longitudinally polarized meson the gluon must also be longitudinal. 
This is, however, forbidden by the gauge invariance as the hadron constituents in 
DA are nearly on-shell. This argument has been known for some time [^. Still, as 



it has a somewhat hand waving flavor (the meaning of being "nearly on-shell" is not 
quite clear), we have formally checked this statement by allowing for longitudinal 
structures in (|T2| , |T3l) and respective extending a set of conformal constraints (see 
Eqs.([T9|- pT]) below). A corresponding system of equations turns out inconsistent. 



which means that the longitudinal polarization in (|T2|,|T3|) is not possible. It is then 
easy to see that the Lorentz structure displayed in Eqs.(p!^,([T3|) is the most general 
decomposition to the 0(twist-5) accuracy. In the early analysis of Ref. ||^ only 
the twist 3 structures $(^),\['(^) were retained on the ground of the fact that they 
correspond to leading power corrections in exclusive amplitudes. They have further 
been analysed within the QCD sum rules technique. We are going to extent the 
analysis of the DA's (|l^,|l^) by a completely different method which is also able to 
fix the "twist 4" 0DA's $(2)^ $(3)^ vl>(2)^ vi>(3)_ Besides of obvious practical interest, we 
find that our new results for the "twist 4" DA's may have some principal importance 
as well. The point is that, as will be shown below, there is a very nontrivial interplay 
between the expansions in twist and conformal spin in the problem at hand. Non- 
leading "twist 4" DA's have lower in comparison to twist 3 DA's conformal spins, 
i.e. are given by poljTiomials of lower dimensions. This observation may imply a 
breakdown of OPE for some processes at the nonleading twist level (see Sect. 5). 
In what follows we restrict our study to the asymptotic DA's, i.e. work to the 
leading order in the conformal spin expansion. A calculation of first pre-asymptotic 
corrections is both possible and desirable. Unfortunately, it is notoriously more 
difficult in comparison to the pion case due to a larger number of possible Lorentz 
structures in matrix elements corresponding to conformal corrections. 

After these lengthy remarks we proceed to the conformal analysis of the DA's 



(|T^),(|T^). As has been explained in the previous Sect. 2, to this end we need to 
construct two light like vectors x^, x^. As for the first one, it is most conveniently 



chosen to be just the light cone coordinate which enters the definitions (12), ( 
(this justifies the identical notations for both). A form of the second light like 
vector bit less obvious. The choice = q^/{qx) which was used in Ref.Q 

for the analysis of the pion DA's clearly does not work in the case at hand as now 
g2 ^ 0. The correct form of the vector for a massive vector meson with arbitrary 
polarization is given by 

= ^ , X = (16) 

supplemented with [e'^^^e^-^^) = 0. It is easy to check that with such choice the 

§Here we put the quotation marks in order to remind that these DA's give rise to further power 
corrections in exclusive amplitudes in comparison to those produced by the DA's One 
has, however, to bear in mind that in the operatorial sense they contain contributions of twist 3, 
as can be seen from Eqs. (0-n9h. 
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projection formulae (]^),(0) and the whole analysis of Sect. 2 remain valid in the 
massive case. 



Selecting different spin components in Eq.([T^) according to (PD,(|9|),([T^) we obtain 
the following system 

(0|(i(-x)7i5(G'*.(fx)M(x)|p(g)) = -m{qx) j D{aiqx)^^^^ (17) 

3 

{0\d{-x)-f,igG4vx)uix)\piq)) = A: / Dia.qx) ($(2) - - $(3) 



m 



{0\d{~x)j^igG^,{vx)u{xMq)) = (e(^)eW)^ ^ D(a,gx)<l>(3) (19) 

{0\di-x)^^tgG±.{vx)uixMq)) = (20) 

plus analogous equations for the set ([13|) which are obtained from (|1^ - ^0]) by 
the substitution igC — > 'j^gG , We have used the fact that for the 

transversely polarized p-meson (e^'^'^x) = for = 0. We will now study separately 
each of the Eqs.(|T7| -p^. Let us start with the first couple of equations. Using the 
general formula (|TT]) we write |^ 

{0\d(—x)'j.igG^.{vx)u{x)\p{q)) = 120ai J D(aiqx)aia2a3{l + . . .) 

(0|d(-x)7.755fG^,.(fx)M(x)|p(g)) = 120bi J D{aiqx) 010203(1 + .. .) (21) 

We introduce here matrix elements of local operators 

{0\dj^igG^,u\piq)) = qaiei^'^Qu - £^^\^.)Al + (el^^^S^, - el^^5^^)m^A2 
{0\d^^^,gG^,u\piq)) = qa{e^^\u-e^^\^)Bi + {e^^^5^,-e^^^5^^)m^B2 (22) 

It is clear that Ai = A2 = Q hy the G-parity, while Bi and B2 are both non-zero 
and will be calculated below. It is then easy to express the normalization constants 
Oi , hi of Eq.(^) in terms of Ai and Bi of (p^). In this way we find 

= 

= 120fiiaia2a3 (23) 

A second couple of constraints is obtained analogously from Eq.(|l^) and its coun- 
terpart for . It reads 

$(1) _ $(2) _^ $(3) ^ 

^(1) _ $(2) + ^(3) = e(5^ + ^2)03 (24) 

The last Eq . (|19D needs more care. Note that the quark spins in (19) are not fixed 
by the transverse gamma matrix. Using the identity 

7± = -^(7-7±7* + 7*7±7-) (25) 



^Here and in what follows ellipses stand for higher conformal spin contributions. 
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we 

line X 



consider two auxiliary DA's with fixed projections of the spins of quarks on the 



{0\d{-x)'y.-f_L'yJgG_L^{vx)u{x)\p{q)) 
{0\d{-x)-f^-f±-f.igG_L^{vx)u{x)\p{q)) 



Qa J D{aiqx)a2{l + ■ ■ \ 
—6a / D{aiqx)ai{l + ■ 



(26) 



The normalization factor a is the same in both expressions by the G-parity. One 
can see that only the antisymmetric combination of gamma matrices contributes in 
Eg. (p6D ■ A little algebra then yields a = B2m^/{qx). After a similar calculation for 
the axial channel we obtain 



$(3) 
^(3) 



3^2 (ttl - ^2) 
3^2(1 -«3) 



(27) 



Using Eqs.( 
totic DA's 



, 0, we finally arrive at the following expressions for the asymp- 



$(1) 

$(2) 
^(3) 





$(3) = 352(ai - aa) 
120-Bia;ia2«3 

3^2(1 - 3^3) + 65ia3(20aia2 - 1) 
3^2(1 - «3) 



(28) 



Note in passing that a peculiarity of restrictions imposed by the conformal invari- 
ance can be seen e.g. in the expression for which does not correspond to any 
simple minded homogenious polynomial in Constants -81,-82 can be calculated 
by the standard QCD sum rules method. In fact, their linear combination Bi +3B2 
(corresponding to the convolution with 6au in Eq.(|2^ )) was found a long time ago 
p!9| . In analogy with the method used in |T^ we find convenient to calculate the 
non-diagonal correlation function 



T{q) = I / dxe^''^(0|T{J(x)7a75^G'^.w(a;)M(0)7AC^(0)}|0) 



(29) 



The matrix element of interest appears as the lowest intermediate hadron state 
contribution to the imaginary part of (^). Matching this expression to a QCD 
calculation at q^ — > —00 we obtain two sum rules for parameters Bi and i?2- The 
Borel transformed sum rules look as follows (a subtraction of continuum from the 
perturbative term is implied): 



Br 
B2 



1 



mpfp 
1 

mpfp 



nip 



24 



M'' + —{:^G')-—7ra 



727r3 ' 24' TT ■ ' 27' 



M2 



+ 



(30) 
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Unfortunately, the first of the sum rules in ( pOj ) is not stable as the condensate terms 
enter there with the same signs. Therefore we use only the second of Eqs. (0) to find 
the number B2 and fix the other parameter Bi from the known value Bi + 3B2 — 
1.85 ■ 10~^ GeV^ which was obtained in Ref.[0 from a similar (stable) sum rule 



written directly for this combination. Numerically this procedure yields 

Bi ~ 0.8 • 10"^ GeV^ , B2 ~ 0.35 ■ 10"^ GeV'^ (31) 

These numbers provide normalizations to the complete system of asymptotic DA's 
(p8|). A comparison between our results (^) and the set of DA's suggested in Ref. 
will be carried out in the next section. Here we would like to end up the discussion of 
the p-meson DA's with a remark that there exists one more series of DA's generated 
by the matrix element 

{0\d{-x)af,^igGaf3{vx)u{x) \p{q)) (32) 



which corresponds to a longitudinally polarized vector meson [20| as there the total 



quark helicity is ±1 and the gluon helicity is =f1. These DA's cannot contribute the 



diffractive amplitude of |n| owing to odd number of gamma matrices in a trace. 
This justifies the claim made in Introduction on the absence of three particle DA's 
in the diffractive electroproduction amplitude of Ref . [H, 



4 Discussions 

Our results ( pH]) should be confronted with the alternative set of DA's suggested in 
Ref.i : 

$(1) = 7!/3^(ai - a2)aia2al , J Da{ai - 02)^^^ = /g^ - 0.25 • 10"^ GeV^ 
= 360jf^aia2al , J Da<i/\p = f^^ ~ 0.6 ■ 10^^ Q^y2 ^33^ 

As the DA's in (|2|) and (g^) are incomparable (any nontrivial is pre- 
asymptotic within a systematic expansion in the conformal spin and thus is beyond 
our asymptotic accuracy), we concentrate on a comparison between the DA's 
First we note that the normalization parameter Bi of the DA's ^^^^ in (^) is consis- 
tent within a typical for the QCD sum rules method accuracy with the normalization 
f^p in (|33|). However, the functional forms of "^^^^ in (pHj) and (|33D are different by an 
extra power of 03. A source of this discrepancy is not quite clear to us. It can be due 
to different definitions of twist in this paper and Ref. 0. Our \E''^^^ corresponds to 
matrix elements of twist 3 operators, while the authors of have in fact singled out 
symmetric parts of operators without subtraction of traces, as just these quantities 
determine power corrections in exclusive amplitudes. Another possible and related 
reason is that when (as in Ref.|Q) the DA is calculated by studying a correlation 
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function of currents with derivatives together with a symmetrization procedure, the 
imaginary part becomes a combination of contributions of different resonances in 
different Lorentz structures. In this case the DA calculated from the asymptotic 
loop may correspond to some "average" hadron (i.e. an admixture of a hypothetic 
exotic meson with J^^ = [|19| is not ruled out). We also note that from the 
viewpoint of the conformal spin expansion the polynomial 0:10203 could arise from 
the next-to-leading term 010203(1 — 803) with j = 7/2, while our asymptotic "^^^^ 
has J = 3 HI]. A coefficient in front of this conformal correction is fixed by the p- 
meson matrix element of the j = 7/2 operator d'j.'j^gD.G^:.u. On the other hand, 
it can be readily seen that this matrix element is exactly zero for the transverse 
polarization of the p-meson. Therefore, it seems unlikely to have a multiplicatively 
renormalizable combination of conformal operators giving rise to the combination 
01O20I and having an anomalous dimension lower than that of the j = 3 operator 
d'j.j^gG^.u which corresponds to the polynomial 01O2O3. To summarize, we see 
a few possible explanations for the difference by a power of 03 between the DA's 
^(1) in (ID and (H), but for its precise identification a further analysis is needed. 
While the choice (^3]) may be ultimately true for the pQCD method of Ref. 0, it 
seems that it is the DA \E'*^^^ of Eq.(28) that has to be used in the light cone QCD 
approach to "soft" exclusive processes 0, since there one deals directly with the 
matrix elements ([I^JT^). 



Now we would like to discuss the non- leading "twist 4" DA's $(2), $(3)^ ^(2)^ ^(3) 
in (p8|) which hitherto have not been considered in the literature. Somewhat surpris- 
ingly, we have found that though these DA's correspond to further power corrections 
in comparison to those produced by the twist 3 DA's their conformal spins 

are lower (j = 2 for $(2), $(3)^ g^^^^ while '^^^^ is a combination of j = 2 and 
j = 3 pieces). Respectively, they are given by polynomials of lower dimensions 
than the asymptotic DA \E''^^) in (^). On the other hand, we remind that typically 
DA's enter exclusive amplitudes through some integrals like / $(oi)/o|, cf. Ref. 0. 
Powers of Oj in the denominator originate from quark and gluon propagators and 
are therefore the same in terms corresponding to three particle DA's of twist 3 and 
4. Moreover, in exclusive amplitudes including three particle DA's of the p-meson 
calculated in |^ corresponding integrals contain logarithmic end-point divergences. 
Our result (^) then implies that if no special cancellations occur (we see no reason 
why they should), in such a situation an attempt to evaluate the twist 4 contribu- 
tion will give rise to power infrared divergences |^ . This would mean a breakdown 
of OPE for this process and invalidation of a calculated twist 3 contribution. This 
observation can be of wider sense and imply that OPE may be inapplicable in exclu- 

II A first possible non-zero term in the DA corresponds to the polynomial (ai — a2)aia2a3 
with j = 7/2. The same power of seems to be missing there as well. 

** Another possibility (V. Braun, private communication) is that power infrared divergences 
may cancel by virtue of equations of motion so that only logarithmic end-point singularities will 
survive. This scenario seems unlikely as the DA's (|8|) are not related by the equations of motion. 
Dynamical calculations are needed to see whether cancellations of infrared divergences take place 
in physical amplitudes. 
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sive reactions beyond the leading twist order. This conjecture does not contradict 
to anything we currently understand about exclusive processes. To illustrate this 
point, we consider two examples. 

1. In the most well studied case of the pion electromagnetic form factor the short 
distance regime breaks down at the 'VjQ'^ order as there the soft Feynman-type 
mechanism p[ brings in the same contribution as a twist 4 part of the hard 
rescattering amplitude (we leave aside the mild Sudakov suppression in the former). 

2. For the 7p7r transition form factor, which is a nonleading twist process related 
to the helicity flipping, the hard rescattering diagram yields the asymptotic l/Q^ 
behavior. However, the soft Feynman contribution gives the same 1/Q^ [0], and 
therefore the short distance regime in this process breaks down already at the lead- 
ing order in l/Q^. If the above scenario were to occur, it would just signal on this 
breakdown, which is otherwise not seen in the leading hard rescattering amplitude 
taken on its own (logarithmic end-point divergences are relatively harmless as they 
can be cured by the Sudakov correction factor). 

A common feature in these two examples is that in the perturbative QCD approach 
the short distance regime is implied within the hard gluon exchange mechanism. 
This assumption certainly does not hold at the nonleading twist level as the 1/Q^ 
parts of the hard amplitudes are just the radiative corrections to the soft Feynman- 
type contributions which have nothing to do with short distances. Thus it is not 
unreasonable to expect divergences at the next order of OPE for these processes as 
a manifestation of the same fact. It would be challenging to see whether there exist 
physical motivations to expect a breakdown of OPE in other and, in particular, 
higher twist exclusive processes like ip ^ pn which seems suggestive according to 
formulae given in Ref. f^. To avoid possible misunderstanding, we should stress that 
in the light cone QCD sum rule approach to exclusive processes the situation is 
different as there the light cone dominance is enforced from outside by an analytic 
continuation to the deep Euclidean region. Therefore, we do not expect a similar 
pattern of a breakdown of OPE in this case, though this possibility cannot be ruled 
out. 
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